Breakdown of a topological phase: Quantum phase transition in a loop gas model with tension 
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We study the stability of topological order against local perturbations by considering the effect of a magnetic 
field on a spin model - the toric code - which is in a topological phase. The model can be mapped onto a 
quantum loop gas where the perturbation introduces a bare loop tension. When the loop tension is small, the 
topological order survives. When it is large, it drives a continuous quantum phase transition into a magnetic 
state. The transition can be understood as the condensation of 'magnetic' vortices, leading to confinement of 
the elementary 'charge' excitations. We also show how the topological order breaks down when the system is 
coupled to an Ohmic heat bath and discuss our results in the context of quantum computation applications. 



Topological phases are among the most remarkable phe- 
nomena in nature. Although the underlying interactions be- 
tween electrons in a solid are not topologically invariant, their 
low-energy properties are. This enhanced symmetry makes 
such phases an attractive platform for quantum computation 
since it isolates the low-energy degrees of freedom from local 
perturbations - a usual cause of errors Tractable theo- 
retical models with topological phases in frustrated magnets 
jjlyJI. Josephson junction arrays ylcl, or cold atoms in traps 
|5[ have been proposed. However, such phases have not, thus 
far, been seen experimentally outside of the quantum Hall 
regime. Is it because topological phases are very rare and 
these models are adiabatically connected only to very small 
regions of the phase diagrams of real experimental systems? 

In this paper, we take a first step towards answering this 
question. We begin with the simplest exactly soluble model 
of a topological phase 1 1], whose Hamiltonian is given below 
in Eq. Q. This model describes a frustrated magnet with four- 
spin interactions similar to cyclic ring exchanges. It is closely 
related to the quantum dimer model 1 6] for frustrated magnets, 
which can be realized in Josephson junction arrays H. We 
consider perturbations of the soluble model that, when suffi- 
ciently large, drive the system out of the topological phase. 
The question is, how large? A small answer would imply that 
this topological phase is delicate and occupies a small por- 
tion of the phase diagram. This might explain the paucity 
of experimentally-observed topological phases. Instead, we 
find that 'sufficiently large' is a magnetic field h of order one 
(h c « 0.6) in units of the basic four-spin plaquette inter- 
action. Our numerical simulations demonstrate, for the first 
time, several key signatures of the phase transition out of the 
topological phase, including the finite-size degeneracy split- 
ting of the topological sectors, the condensation of 'magnetic' 
excitations, and the confinement of 'electric' charges. 

We also consider perturbing the system by coupling it to an 
Ohmic heat bath. When coupled to such a bath, a quantum 
mechanical degree of freedom can undergo a transition from 
coherent to incoherent behavior |7]. Recently, the effects of 
such a coupling on quantum phase transitions, at which di- 



vergent numbers of quantum mechanical degrees of freedom 
interact, have also been studied 1 8]. In both contexts, the cou- 
pling to the heat bath tends to make the system more classical. 
Coherent quantum oscillations are suppressed, while broken 
symmetry phases - which are essentially classical - are stabi- 
lized. A topological phase is quantum mechanical in nature. 
We find that coupling the heat bath to the kinetic energy, i.e. 
plaquette flip operator, does not destroy such a phase. How- 
ever, when the dissipation is strong, the gap becomes very 
small, and the topological phase may be too delicate to ob- 
serve or use at reasonable temperatures. On the other hand, if 
the heat bath is coupled to the classical state of each plaque- 
tte, the topological phase is destroyed through a Kosterlitz- 
Thouless transition at a dissipation strength of order one. 

These results can be recast in quantum information lan- 
guage: the ground states in different topological sectors are 
the different basis states of an encoded quantum memory. 
Quasiparticle excitations are states outside of the code sub- 
space. The stability of the topological phase, as measured by 
an energy gap A within a topological sector (essentially the 
energy cost for a pair of quasiparticles), translates into an er- 
ror rate for topological qubits. At zero temperature, errors are 
due to the virtual excitation of pairs of quasiparticles, assum- 
ing that the system is shielded from perturbations at frequen- 
cies higher than A. Such virtual processes lead to a splitting 
between topological sectors 5E ~ e~ AL / v , where L is the 
system size and v is a characteristic velocity. As the tempera- 
ture is increased, the thermal excitation of particles eventually 
becomes more important and the error rate is ~ e~ f3A l25ll . 
(The actual concentration of excitations leading to unrecover- 
able errors was studied in IMm .) 

The model - We start with the toric code Hamiltonian 01 

h tc =-ay: n n 

v j£vertex(f) V J (Eplaquette(p) 

where the cr, are S = 1/2 quantum spins on the 27V edges 
of a square lattice with N vertices on a torus. Since all terms 
in Hamiltonian Q commute with each other, the model can 
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be solved exactly 1 1]. The ground-state manifold can be de- 
scribed as a quantum loop gas where the loops consist of 
chains of up-pointing a z -spins and the loop fugacity is d = 1. 
On the torus there are four degenerate ground states which can 
be classified by a winding number parity P y / X = Yiiec / °f 
along a cut c x / y in the x— or y— direction. 

Here we study the effect of perturbing the Hamiltonian Q 
with a loop tension which can be introduced either by a lon- 
gitudinal magnetic field or local Ising interaction of the form 

H=H TC -hj2 at - jY > w 

where h (J) is the strength of the magnetic field (Ising inter- 
action). These are the dominant perturbations expected in a 
physical implementation, e.g. in a Josephson junction imple- 
mentation yl 0] they arise from electric potential perturba- 
tions or Coulomb interactions between neighboring quantum 
dots. We discuss this model in the limit of a large charge gap, 
i.e. A 3> B,h, J, where it becomes equivalent to the 'even' 
Ising gauge theory II ill . The low-energy sector has no free 
charges and any state is described by a collection of loops that 
can be obtained from a reference state (e.g. all erf = 1/2) by 
a sequence of plaquette flips. Let us introduce a new plaquette 
spin operator [i p with eigenvalues fj,* = (— l)" p /2 where n p 
is the number of times a given plaquette p has been flipped, 
counting from the reference state. Then of = 2fi z fi z , where 
p and q are the plaquettes separated by the edge i. The pla- 
quette flip term in Eq. Q becomes — 45 ^ p fi p . In the new 
variables, Hamiltonian (|2} becomes equivalent to the trans- 
verse field Ising model (with both nearest and next-nearest 
neighbor Ising interactions) in a basis restricted to loop states. 
Independent of the choice of basis states the system orders 
at a critical field strength (h/B) c = 0.65695(2) determined 
by continuous-time quantum Monte Carlo simulations [12]. 
The transverse field Ising model for the plaquette spins can be 
mapped to a classical (2+l)-dimensional Ising model: 

Wei = -K T Y S p (t)S p (t + At)-K Y, Sp(T)S q (T), 

T,p T,{p,q) 

(3) 

where the imaginary time r is discretized in steps of At and 
S p = 2/ip = ±1. The magnetic field h and the local Ising 




FIG. 1: Left: Illustration of the toric code Hamiltonian Q. Right: 
The elementary excitations above the loop gas ground state are pairs 
of 'magnetic vortices' on plaquettes connected by a string of a z - 
operators (solid line) and 'electric charges' on the vertices connected 
by a string of a x -operators (dashed line). 
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loop tension h 

FIG. 2: (Color online) Magnetization M = (£\ of }/2N versus 
loop tension (magnetic field). The topological phase survives for 
small loop tension where an almost constant susceptibility (see in- 
set) leads to a linear increase of M. Above the critical loop tension 
(dashed line) the system approaches the fully polarized state. 

interaction J introduce a nearest and next-nearest neighbor 
interaction between the classical spins S p (t), respectively. 
Since the exact nature of this Ising interaction does not play a 
role in the following we do not discuss the case of next-nearest 
neighbor couplings in detail. The coupling along real-space 
coordinates is then given by K — ^Ar • h and along imagi- 
nary time by K T = — ^ In [tanh(Ar ■ B)\. The model de- 
scribes the well-known continuous magnetic phase transition 
of the 3D Ising model. For isotropic interactions, K = K T , 
the critical coupling has been determined with high precision 
to be K c = 0.2216595(26) Q. Setting B = 1 this gives a 
critical loop tension h c = 0.58224 with isotropic couplings at 
At = 0.761403. 

The magnetic susceptibility diverges at the transition and 




loop tension h 

FIG 3: (Color online) Excitation gap for magnetic vortex excitations 
(star symbols) versus the loop tension. At the critical loop tension 
the gap closes and the magnetic vortices condense. Right ordinate: 
Rate exp(— /3A) for tunneling processes between topological sectors 
(filled circles) for /3 = 10. 
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FIG. 4: (Color online) Splitting of the energy degeneracy for the four 
topological sectors defined by even or odd winding number parities 
along the x- and y-directions. The intermediate {even-odd / odd- 
even} -parity sector was taken as a reference (dashed line). 



FIG. 5: (Color online) Finite-size scaling of the energy splitting 
between the topological sectors around the quantum phase transition 
for j3 — 10L. At the critical loop tension (crosses) we find power- 
law scaling 5E{L) oc L 2_z with exponent z = 1.42 ± 0.02. 



the magnetization Q^. of )/27V shows a corresponding kink, 
as shown in Fig. [2] Although this is not a symmetry -breaking 
transition, the analogous transition driven by next-nearest in- 
teraction J is a continuous quantum phase transition from a 
topologically ordered quantum state to a broken symmetry 
state 1 14]. The magnetic transition can also be understood 
in terms of the condensation of 'magnetic vortices', plaque- 
ttes with Y[j &j = — !■ While for the original Hamiltonian 
([0, the gap to these excitations is A = 2B, they become 
gapless and condense at the critical loop tension, as shown in 
Fig. [3] The gap has been estimated from measurements of 
the imaginary time correlation length £ T as A oc l/£ r which 
we have calculated applying continuous-time quantum Monte 
Carlo simulations using the ALPS looper code OS)- The 
thermal excitation of pairs of magnetic vortices occurs with 
probability ~ cxp(— /3A), also shown in Fig. [3] 

Topological order - The breakdown of topological order 
at the phase transition can be seen from the energy splitting SE 
between the ground-states for the various topological sectors. 
When winding parities are used as basis states for a quantum 
memory, this splitting causes phase errors. (The absence of 
'electric' charges precludes any transitions between different 
winding parities so bit flip errors cannot occur.) In the topo- 
logical phase, the virtual excitation of quasiparticles leads to 
a small splitting SE oc exp(— AL/v) between the topolog- 
ical sectors. In the classically-ordered phase, on the other 
hand, the energy splitting should scale with L, which cor- 
responds to the energy cost of a loop in the ordered ground 
state. As the winding parity is conserved by imaginary time 
spin-flip operations, we can simulate the system in one of 
the topological sectors by defining an initial spin configura- 
tion that corresponds to the respective limit for large loop ten- 
sion. Fig.|4]shows the calculated splitting for various system 
sizes in the vicinity of the critical loop tension. At the phase 
transition, the behavior qualitatively changes from power-law 
scaling for strong loop tension to an exponential suppres- 



sion in the topological phase for small loop tension as dis- 
cussed above. A more quantitative picture arises from the 
finite-size scaling analysis of the energy splitting 5E(L) be- 
tween the {even-odd}- and {even-even}-parity sectors shown 
in Fig. [5] For the critical loop tension we find a power-law 
scaling SE(L) oc L 2 ~ z with an exponent z = 1.42±0.02. Be- 
low the critical value the scaling turns into exponential scaling 
as expected for the topological phase. 

Confinement transition - For the loop gas Hamiltonian 
([0 the elementary electric charge excitations (end-points of an 
open loop) are deconfined. For strong loop tension, however 
these excitations are expected to become confined, thereby 
eliminating all open loops. We can study this confinement 
transition in our simulations of model by artificially in- 
troducing pairs of electric charge excitations for the sam- 
pled loop configurations. This allows us to measure the con- 
finement length £ c as the square root of the average second 
moment of the distance between the two excitations, which 
for a torus with even extent L has to be normalized by a 
factor 6/(L 2 + 2). The measured confinement length £ c 
shown in Fig. |6] clearly demonstrates that electric charges re- 
main deconfined for the full extent of the topological phase 
and that the confinement transitions occurs simultaneously 
with the magnetic phase transition. At the critical loop ten- 
sion the confinement lengths £ C (L) for various system sizes 
cross which demonstrates that the confinement length di- 
verges with the same critical exponents as the magnetic cor- 
relation length £ and there is only one length scale describ- 
ing the phase transition. This can be seen from the follow- 
ing finite-size scaling argument: We may write the finite-size 
scaling behavior for the confinement length as £ c (r, L) oc 
t~ v j (Lt x ) where v (x) is the critical exponent of the cor- 
relation (confinement) length respectively and r = h — h c . 
With /(L 1/x t) = (i 1 /x T )-"/( J L r x) we obtain for the con- 
finement length £ c (t, L) oc L u / x f(L x / x T). At the phase 
transition the existence of a crossing point £ c (0, L)/L oc 
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FIG. 6: (Color online) Confinement length of two electric charge 
excitations versus loop tension. The confinement transition occurs at 
the same critical loop tension (dashed line) as the magnetic transition. 



FIG. 7: (Color online) Gap and error probability versus dissipa- 
tion strength estimated from the correlation time £ T of the dissipative 
Ising chain J4j. 



[jV/x- 1 = const then implies \ = v. For our model without 
dynamical electric charges, this measure of the confinement 
of test charges is closely related to the calculation of a Wilson 
loop expectation value. In the presence of dynamical elec- 
tric charges, this becomes trickier; Polyakov loops have been 
used as an order parameter for the finite temperature transition 
of the 3D Ising gauge model 1 171. 

Dissipation - Finally, we discuss the effect of dissipa- 
tion when Hamiltonian Q is coupled to an Ohmic heat bath. 
Since our model excludes dynamical electric charges, we do 
not consider coupling a heat bath directly to erf. Instead, we 
first examine coupling a heat bath to so that a 'phonon' is 
created when a plaquette flips. This type of dissipation could 
occur in a Josephson junction model 1 18[ or in a spin model 
through the spin-phonon coupling. The standard procedure 
lll9ll for a linear spectral density ('Ohmic' dissipation) results 
in an effective action for independent Ising chains with long- 
range couplings in an external longitudinal magnetic field 
(which here means parallel to fi x ). We note that as a con- 
sequence of the Lee- Yang theorem 1 20], there can be no sin- 
gularities of the respective partition function at any real non- 
zero longitudinal field, ruling out the existence of a quantum 
phase transition for this model. In particular, this implies that 
the magnetic gap remains finite for any dissipation strength! 

An entirely different behavior arises if dissipation is cou- 
pled such that it stabilizes the 'classical' state of the system. 
Coupling the bath to either of or fjJt stabilizes the classical 
state of a single spin or a plaquette, respectively. We consider 
the latter as it should be more effective at damping quantum 
fluctuations, although we expect the former to have similar 
physics. The same procedure as above then leads to a model 



for decoupled Ising chains given by 

Ha = -K T S p {t)S p {t + At) 

_ a (JL.Y ,S 'p( t )' 5 p( t ') (a\ 

2 ^ \N T ) s in 2 (^|T~r'|) ' K ) 

t<t' ,p v N t I IJ 

where the parameter a measures the dissipation strength. This 
model has been well studied |21] and is known to exhibit a 
Thouless-type phase transition into a classically disordered, 
fluctuationless phase. The critical value a c of this transition 
depends weakly on the cutoff in the long-range interaction; in 
our simulations a c « 0.7. At the transition the magnetic gap 
vanishes, in sharp contrast to the previous case. 

Due to the long range interactions introduced by the dissi- 
pative coupling, the spin-spin correlations will asymptotically 
decay as 1/r 2 12^1 . It is therefore a non-trivial task to de- 
fine a correlation-time and hence to estimate the excitation 
gap. For a < 0.1 one observes an exponential decay of the 
correlation function onto the asymptotic 1/r 2 behavior, and 
subtracting the l/r 2 -contribution allows one to estimate £ r , 
which is found to grow linearly with a. For a > 0.1, this 
procedure can no longer be used. In this region we estimate 
£ r from the asymptotic decay of the correlations, proportional 
to {£, t /t) 2 . These £ T s grow approximately exponentially in 
the region 0.1 < a < a c . Alternatively, one could define a 
correlation-time from the crossover scale where the short-time 
behavior crosses over to the asymptotic 1 /r 2 form. This defi- 
nition in the spirit of a 'Josephson length' I23I1 yields (up to a 
normalization factor) the same results as the £ T extracted from 
the asymptotic decay. The gap estimated from the inverse cor- 
relation function, as well as the error probability is plotted in 
Fig. We find that the error probability remains negligibly 
small below the crossover value a 1=1 0.1. 

Outlook - We have shown that the topological phase 
which governs the toric code model (U actually exists in an 
extended region of phase space around the soluble point. It 
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is stable against deviations of the system Hamiltonian from 
the ideal one and also the coupling of the system to its en- 
vironment. In general, this demonstrates that a system does 
not necessarily have to be particularly fine-tuned to reach a 
topological phase. The paucity of their experimental obser- 
vations may thus be due not to some intrinsic delicateness of 
such phases, but rather to the experimental subtlety involved 
in identifying them. In future work, these conclusions need 
to be tested in other, more exotic topological phases which 
support universal topological quantum computation 12411 . 
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